Z 2 topological term, the global anomaly, and the two-dimensional symplectic 
symmetry class of Anderson localization 
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We discuss, for a two-dimensional Dirac Hamiltonian with random scalar potential, the presence 
of a Z 2 topological term in the non-linear sigma model encoding the physics of Anderson localization 
in the symplectic symmetry class. The Z 2 topological term realizes the sign of the Pfafhan of a family 
of Dirac operators. We compute the corresponding global anomaly, i.e., the change in the sign of 
the Pfaffian by studying a spectral flow numerically. This Z 2 topological effect can be relevant to 
graphene when the impurity potential is long-ranged and, also, to the two-dimensional boundaries 
of a three-dimensional lattice model of Z 2 topological insulators in the symplectic symmetry class. 
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There is strong theoretical and numerical supporting 
evidence for the premise that metal-insulator transitions 
in the problem of Anderson localization with short-range 
correlated disorder can be classified in terms of a few fun- 
damental properties of microscopic Hamiltonians. When- 
ever the motion of the relevant (fermionic) quasiparti- 
cles is diffusive, a non-linear sigma model (NLerM) can 
be derived for each symmetry class. A NLerM encodes 
fluctuations of Nambu-Goldstone bosons (diffuson and 
Cooperons) that are defined on a curved manifold (target 
manifold) determined by the symmetries of a microscopic 
disordered Hamiltonian. 

When the fluctuations are small, the dynamics of the 
Nambu-Goldstone bosons are solely determined by the 
local data of the target manifold such as the metric and 
curvature. However, the global topology of the target 
manifold does have important effects. For example, a 
weak magnetic field relative to the disorder strength leads 
to localization of all states in two-dimensional space. On 
the other hand, a strong magnetic field leads to the in- 
teger quantum Hall plateau transition. This difference 
is captured by the absence or presence, respectively, of 
a topological term in the NLerM corresponding to the 
two-dimensional unitary symmetry class Similarly, a 
random single-particle Hamiltonian that preserves time- 
reversal symmetry but breaks spin-rotation symmetry 
defines the symplectic symmetry class and can be as- 
sociated to a NLerM. In two-dimensional space, Fendley 
Q made the observation that such a NLerM admits a 
topological term. The configuration space of the relevant 
NLerM has two topologically distinct sectors. Configura- 
tions from different topological sectors can thus be given 
different Boltzmann weights in the presence of the topo- 
logical term, thereby realizing the outcome of a metal- 
insulator transition belonging to a universality class dif- 
ferent from the one without the topological term Q . 

We show in this letter that a topological term, encod- 
ing a global (i.e., nonperturbative) anomaly, is realized 
by the Pfaffian of Majorana fermions that originate from 



the problem of Anderson localization defined by the two- 
dimensional two-component Dirac Hamiltonian 

d 



H= er -p + cr o y(r) 



P:=VT, (1) 
idr 



subject to a white-noise and Gaussian correlated random 
scalar potential V(r), V(r) = 0, V(r)V(r') = gS(r—r'). 
Here, r € R 2 , g > measures the disorder strength, and 
<J X y z denote the standard Pauli matrices with cr the 
corresponding unit matrix. In the recent independent 
derivation by Ostrovsky et al. 0, the topological term 
is not presented in the form of the sign of a Pfaffian 
which, however, is essential for it to be interpreted as 
a global anomaly. The fabrication of graphene 0] has 
triggered a renewed theoretical interest in the properties 
of the random Dirac Hamiltonian ^ . We will argue that 
the random Dirac Hamiltonian |T]) can be realized at the 
two-dimensional boundaries and in the low-energy limit 
of three-dimensional lattice models. 

As observed by Ludwig et al. [B| the symmetry 
iUyTL* (—i(Ty) = 7i, which we shall abusively call time- 
reversal symmetry (TRS) as it might not necessarily 
realize the TRS of the underlying microscopic model, 
puts the random Dirac Hamiltonian |T|) in the two- 
dimensional symplectic symmetry class. The disorder in 
the Dirac Hamiltonian |T]) is thus expected to yield weak 
anti- localization corrections to the mean conductance [f| . 

The link between single-particle random Hamiltoni- 
ans and NLerM comes about when setting up a gener- 
ating function for the mean value taken by the product 
of the Green's functions. Using the standard machin- 
ery of the replica trick [3], the disorder-averaged prod- 
uct of the retarded and advanced Green's functions of 
Eq. JT]) can be obtained from the Grassmann path inte- 
gral Z := J D yj, ip\ exp (— J cPr C) , where 

£ = 4>a (ff -P5ab + "7 A ab) i'b - n^a^aMb, ( 2 ) 



and with the Majorana condition ip a := 
Grassmann integration variables (a = 1 , 
mation over repeated indices with a, b 



■0a «cr y , on the 
..,4iV r ). Sum- 
= 1,...,4AL is 
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implied. For any a = 1, . . . , 47V r , the Pauli matrices act 
on the two components of the spinor ip a . The 47V r x 4iV r 
diagonal matrix A = diag (l2Ar r , — l2JV r ) distinguishes two 
sectors coming from the retarded and advanced Green's 
functions. It is multiplied by the infinitesimal positive 
number ?y. The replica limit 7V r — > is understood if the 
effective partition function is purely fermionic, while the 
choice N r = 1 is appropriate when ((2]) is regarded as the 
fermionic sector of a supersymmetric representation of 
the product of the retarded and advanced Green's func- 
tions. The 4-fermion interaction in @ originates from 
the integration over the static random field V(r). 

When r] = 0, the replicated action is invariant under 
a 0(4N T ) transformation acting on the (time-reversal) ® 
(retarded/advanced) ® (replica) indices. The small imag- 
inary part of the energy rj lifts the degeneracy between 
the retarded and advanced sectors, thereby reducing the 
0(47V r ) symmetry down to 0(27V r ) x 0(2N r ). The 
Nambu-Goldstone modes associated with this symmetry 
breaking are described by a bosonic matrix field Q(r) liv- 
ing in the cosct space G/H = 0(4JV r )/0(2iV r ) x 0(2N r ). 
They emerge after performing the Hubbard-Stratonovich 
decoupling of the interaction by a 0(4N 1 .) matrix field 
and then freezing its massive modes (i.e., all modes other 
than the Nambu-Goldstone modes). Their interactions 
are governed by the effective local field theory 



Z eS = V 



[^]V[Q] e -/^[A+^tr(gg-)], 



(3a) 



C { := ip D[Q] tp = ip T D[Q] ijj. 
Here, we have introduced the Majorana kernel 

(5[Q]) Qb :=io y (a -p5 ab - iAQ ab + i V A ab ) 

=:ia y (D[Q]) ab , a, b = 1, . . . , W T . 



(3b) 



It is skew symmetric, (icr y D[Q\j = —ia y D[Q], in the 
background of Q S G/H, i.e., 



Q 2 =h N , Q = <3 T ! trQ = 0. 



(3c) 



The real number A, which represents the imaginary 
part of the self-energy caused by disorder, is deter- 
mined by solving the self-consistent equation ir/g = 
In [l + (a A)™ 2 ] where a is the short-distance cutoff. 

For smooth spatial variations of the Q-field starting 
from A, the spectrum of D[Q] has a gap of order A 
at the band center. We can then ignore rj and in- 
tegrate the Majorana fermions out. An effective ac- 
tion for the Nambu-Goldstone modes follows. This in- 
tegration, when properly regularized, defines the Pfaf- 

fian / V V] exp (-/ d 2 r £ f ) = Pf (£>[Q]^j . A possi- 
ble regularization that preserves TRS starts with wrap- 
ping momentum space around the two-torus T 2 . Mo- 
mentum is then quantized, k = 2itn u /L where fi = 1,2, 



with the level spacing 2tt/L controlled by the long- 
distance cutoff L. Momentum is made countably finite, 
= —(N - l)/2,-(N- l)/2+ l,...,(iV- l)/2 where 
fx = 1,2 with the help of the short-distance cutoff a . 
The ratio N = L /a of the ultraviolet to infrared cutoff 
is chosen an odd integer for convenience. To complete 
the definition of the regularized Pfaffian we note that its 
square is the determinant DetD[Q]. This determinant 
can be rewritten 

Det D[Q] = Dct (ia y D'[Q]) = Dct D'[Q] (4a) 

with the Dirac kernel (a,b = 1, . . . , 4A^ r ) 

(D'[Q]) ab = <r.p6 ab + <T z AQ ab . (4b) 

The real- valued spectrum of D'[Q] comes in pairs of non- 
vanishing eigenvalues {— A^,+A^} labeled by the index i 
running over some countably finite set. For some refer- 
ence configuration Q, we define 



ViD[Q] = \{K 



(5) 



where A, S R represents either one of — A^ or +A^ for each 
i. It follows that argPf D[Q] £ {0,tt}. The sign of the 
Pfaffian is protected by the spectral gap oc A under any 
infinitesimal change of Q. This protection does not ex- 
tend to all configurations. Indeed, the second homotopy 
group of the target manifold is non-trivial, 



So 



n 2 (G/H) 



for all N x > 1, 



ZxZ, for 7V r = 1, 



(6) 



and, as we are going to show explicitly, there is a phase 
difference of ir between the two Pfaffians (O when evalu- 
ated at two Q' s belonging to distinct Z 2 topological sec- 
tors. This ambiguity in fixing the sign of the Pfaffian 
over all configurations, starting from some reference one, 
is reminiscent of the SU(2) global anomaly in four space- 
time dimensions that follows from 7r 4 (SU(2)) = Z 2 [8[. 

We turn to the construction of two families of Q-ficlds 
that differ in the sign of the Pfaffian ((5]). To this end, we 
first compactify two-dimensional Euclidean space r G R 2 
by wrapping it once around the two-sphere S 2 . On the 
two-sphere parametrized by the polar — ir/2 < 9 < +tt/2 
and azimuthal < <f> < 2tt angles, we define, following 
Weinberg et al. [J the family of fields 



Q 



l 2N _ 2 
q k (6,<t>) 
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2JV, -2 



?*(M) = 

RM = 



cos ( 



sin 6Rl(cj)) 

cos k(f> sin k<f> 
- sin kcf> cos kcf> 



sm6R k (<j>) 



(7) 
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(a) (Q t ,Q f ) = (A, 0*=,,) 

1.0 



(b) (Qi,Q f ) = {Q h=0 ,Q^) 
1.0 




(c) (Q„Q f ) = (Qk= ,Qk= 2 ) (d) (Qi,Qf) = (A,Q fc=1 ) 

1.0 




FIG. 1: The energy eigenvalue spectrum in the vicinity of the 
band center for the Dirac kernel D'[Q t ], Eq. ()4b|) , is computed 
numerically as a function of the parameter < t < 1 for 
(aoA) -1 = 1 and N = 11. The field Q t interpolates between 
Qi when t — and Qf when t — 1. (e) and (f): Same as 
in (c) and (d), respectively, except for the addition of a small 
perturbation so as to lift any accidental or quasi degeneracies. 



labeled by the index k G Z. One verifies 
that the Chern integer (Pruisken term) Ch[Q fe ] := 
(16ni) 1 J S 2 d?r e M „tr \Qk^^Qk^vQk\ vanishes for any 
k G Z. This is expected as TRS holds. On the other 
hand, we are going to argue that 

sgnPf (5[Q fe ]) = -sgnPf (5[Q fe+1 j) , k G Z. (8) 

Any clement Q G G/i? is also an element of the larger 
coset space U(4A r ) /U(2iV r ) x U(2A r ) . From this point of 
view, the Chern integer is the signature of the second ho- 
motopy group tt 2 (U(4JV r )/U(2JV r ) x U(2iV r )) = Z. As 
the Chern integer is also the phase of the fermion de- 
terminant of the Dirac kernel £>'[[/] (through the chiral 
anomaly) in the background of U G U(4JV r )/U(2JV r ) x 
U(2A r ), the fact that the Chern integer vanishes is a 
consequence of the positivity of the determinant (|4a|) . 
While the Chern integer is blind to the second homotopy 
group ((6J, the Pfaffian is sensitive to it. 

We give a "numerical proof" of Eq. ([8|) in Fig. Q] by 
showing 4 evolutions of eigenvalues of the Dirac ker- 
nel jib)) evaluated at Q t := (l — t)Q i +tQj : as a function 
of < t < 1. Here, the initial, Q i; and final, Qf, field 
configurations belong to G/H [8[. Although Q t is not a 
member of G/H for < t < 1, it remains real- valued, 
symmetric, and traceless. Consequently, the spectrum of 
D'[Q t ] is symmetric about the band center. Configura- 



tions Q i and Q ^ have Pfaffians of opposite signs whenever 
an odd number of level crossing occurs at the band cen- 
ter ( "spectral flow" ) during the devolution of the Dirac 
kernel D'[Q t ]. This is accompanied by the closing of 
the spectral gap of D'[Q t ] by an odd number of pairs 
(— A-,+A-) as t interpolates between and 1. The spec- 
tral i-evolution is obtained numerically using the regu- 
larization of the Dirac kernel D'[Q t ] in momentum space 
as described above Eq. (j4]) [T3|. We show with Fig. QJa) 
that A (the uniform configuration) and Q k= Q belong to 
the same Z 2 topological sector as the spectral gap never 
closes under the evolution of the spectrum: Pf (D[A}) and 
~P{(D[Q k=0 \) share the same sign. We show with Fig. 
[TJb) that Q k=0 and Q k=1 belong to different Z 2 topolog- 
ical sector as level crossing at the band center takes place 
for a single pair of levels: Pf(D[Q k=0 \) and Pf(D[Q k=1 ]) 
differ by their sign. We show with Figs.[ljc) and (e) that 
Q k=Q and Q k=2 belong to the same Z 2 topological sec- 
tor as level crossing at the band center takes place for 2 
pairs of levels [3]: Pf(-D[Q fc=0 ]) and Pf(£>[Q fc=2 ]) share 
the same sign. Finally, we show with Fig. [Hd) and (f) 
that A and Q k=1 belong to different topological sectors 
as level crossing at the band center takes place for 3 pairs 
of levels [IH: Pf(5[A]) and Pf(5[Q fe=1 ]) have opposite 
signs. 

From now on, we shall assign the topological quantum 
number (1) to all configurations belonging to the same 
homotopy class as Q k with k even (odd). If so, the effec- 
tive action (|3aj) can be approximated by [T2I ] 



^topolo 

NLctM — 



V\Q\ (-1)»W1 e~ s ^ 



(9) 



where S[Q] is the usual local action for the NLcM on 
G/H and n[Q] = 0,1 is the Z 2 topological quantum 
number of Q. The topological term has its origin in 
the Pfaffian arising from Majorana fermions, i.e., the 
global anomaly. When the bare NLerM coupling con- 
stant is small (when the bare conductivity is large), the 
effect of n[Q] = 0, 1 on the renormalization group flow is 
small. Weak antilocalization, as it occurs in the absence 
of the topological term, is then expected. The effect of 
n[Q] = 0, 1 is more important in the strong coupling 
regime (when the bare conductivity is small) 0. Nu- 
merical simulations HI, 14 1 of the random Dirac Hamil- 
tonian ((T|) suggest that the topological term makes the 
two-dimensional symplcctic stable metallic fixed point 
the only fixed point. 

As is well known from lattice gauge theory [jjj . the 
fermion doubling problem forbids the emergence of an 
odd number of massless spinors from two-dimensional 
tight-binding models with conserved electric charge, 
TRS, locality, and the translation symmetry of a reg- 
ular lattice. That is, the NLcrM ((9]) cannot emerge from 
microscopic electronic models on the square or honey- 
comb lattices [l6|. However, as pointed out by Ando and 



4 



Suzuura FJj, graphene approximately realizes the sym- 
plectic symmetry if the potential range of static point 
impurities is much larger than the lattice spacing. 

One way out of the fermion doubling problem is 
to consider d-dimensional space as the boundary of 
(d + l)-dimensional space 0, [lj)[ and to study the 
physics of localization at this boundary. In this con- 
text, the effect of disorder on the two-dimensional bound- 
ary states of three-dimensional Z 2 topological insulators 



2 ll . l22j could yield microscopic realizations of the 



two-dimensional symplcctic class with a Z 2 topological 
term. The situation is here similar to that in a quasi-one- 
dimcnsional tight-binding model belonging to the sym- 
plcctic symmetry class. When the number of Kramers 
doublets propagating in the wire is even, all states are ex- 
ponentially localized j23[. When it is odd, one Kramers 
doublet remains extended 2J, |25|. The former case is 



always realized in a quasi-one-dimensional tight-binding 
model [23[ , whereas the latter case can be realized at the 
edges of two-dimensional Z 2 topological insulators [26| . 
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